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Acoustic metafluids are defined as the class of fluids that allow one domain of fluid to acoustically 
mimic another, as exemplified by acoustic cloaks. It is shown that the most general class of acoustic 
metafluids are materials with anisotropic inertia and the elastic properties of what are known as 
pentamode materials. The derivation uses the notion of finite deformation to define the transforma- 
tion of one region to another. The main result is found by considering energy density in the original 
and transformed regions. Properties of acoustic metafiuids are discussed, and general conditions 
are found which ensure that the mapped fiuid has isotropic inertia, which potentially opens up the 
possibility of achieving broadband cloaking. 



PACS numbers: 43.20.Fn, 43.20Tb, 43.40Sk, 43.35. Bf 



I. INTRODUCTION 

Ideal acoustic stealth is provided by the acoustic cloak, 
a shell of material that surrounds the object to be 
rendered acoustically "invisible" . Stealth can also be 
achieved by "hiding under the carpet" ^ as shown in 
Fig. 1. A simpler situation but one that displays the 
essence of the acoustic stealth problem is depicted in 
Fig. 2. The common issue is how to make one region 
of fluid acoustically mimic another region of fluid. The 
fluids are different as are the domains they occupy; in 
fact the mimicking region is typically smaller in size, it 
can be viewed as a compacted version of the original. 

The subject of this paper is not acoustic cloaks, or 
carpets, or ways to hide things, but rather the type of 
material necessary to achieve stealth. We define these 
materials as acoustic metafluids, which as we will see can 
be considered fluids with microstructure and properties 
outside those found in nature. The objective is to derive 
the general class of acoustic metafluids, and in the pro- 
cess show that there is a closed set which can be mapped 
from one to another. Acoustic metafluids arc dcflned as 
the class of fluids that (a) acoustically mimic another re- 
gion as in the examples of Figs. 1 and 2, and (b) can 
themselves be mimicked by another acoustic metafluid 
in the same sense. The requirement (b) is important, 
implying that there is a closed set of acoustic metafluids. 
The set includes as a special case the "normal" acous- 
tic fluid of uniform density and bulk modulus. Acoustic 
metafluids can therefore be used to create stealth devices 
in a normal fluid. But, in addition, acoustic metafluids 
can provide stealth in any type of acoustic metafluid. 
The reciprocal nature of these fluids make them a natu- 
ral generalization of normal acoustic fluids. 

The acoustic cloaks that have been investigated to date 
fall into two categories in terms of the type of acoustic 
metafluid proposed as cloaking material. Most studies, 
e.g.^""^, consider the cloak to comprise fluid with the nor- 
mal stress-strain relation but anisotropic inertia, what 
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FIG. 1. Three ways to acoustically hide something: envelope 
it with a cloak (top); hide it under a carpet (middle); or hide 
the object in the corner. In each case the acoustic metafluid 
(red) emulates the acoustic properties of the hidden region as 
if it were filled with the exterior fluid (blue). 



we call inertial cloaking. Particular realization of incrtial 
cloaks are in principle feasible using layers of isotropic 
"normal" fluid^^"^^; the layers are introduced in order 
to achieve an homogenized medium that approximates a 
fluid with anisotropic inertia. An alternative and more 
general approach^' ^-'^ is to consider anisotropic inertia 
combined with anisotropic elasticity. The latter is intro- 
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FIG. 2. The top shows waves in a semi-infinite medium s > 0. 
The wave incident from the right reflects from a perfectly 
reflecting boundary at a; = 0. The lower figure shows the 
same medium in a: > & with the region < a; < 6 replaced by 
a shorter region of acoustic metafluid. Its properties are such 
that it produces a perfect acoustic illusion or "mirage". 



duced by generalizing the stress strain relation to include 
what are known as pentamode^ elastic materials®'^"'^^. 
Clearly, the question of how to design and fabricate 
acoustic nietafluids remains open. The focus of this paper 
is to first characterize the acoustic metafluids as a general 
type of material. In fact, as will be shown, this class of 
materials contains broad degrees of freedom, which can 
significantly aid in future design studies. 

The paper is organized as follows. The concept of 
acoustic metafluids is introduced in §11 through two 
"acoustic mirage" examples. The methods used to find 
the acoustic metafluid in these examples are simple but 
not easily generalized. An alternative and far more pow- 
erful approach is discussed in §111: the transformation 
method. This is based on using the change of vari- 
ables between the coordinates of the two regions com- 
bined with differential relations to identify the metafluid 
properties of the transformed domain. Leonhardt and 
Philbin^^ provide an instructive review of the transfor- 
mation method in the context of optics. The transforma- 
tion method docs not however define the range of mate- 
rial properties capable of being transformed. This is the 
central objective of the paper and it is resolved in §IV by 
considering the energy density in the original and trans- 
formed domains. Physical properties of acoustic metaflu- 
ids are discussed in §V, including the unusual property 
that the top surface is not horizontal when at rest un- 
der gravity. The subset of acoustic metafluids that have 
isotropic inertia is considered in §VI, and a concluding 
summary is presented in §VII. 



II. ACOUSTIC MIRAGES AND SIMPLE METAFLUIDS 

The defining property of an acoustic metafluid is its 
ability to mimic another acoustic fluid that occupies a 
different domain. The simplest type of acoustic illusion is 
what may be called an acoustic mirage where an observer 
hears, for example, a reflection from a distant wall, but in 
reality the echo originates from a closer boundary. Two 
examples of acoustic mirages are discussed next. 



A. ID mirage 

Consider perhaps the simplest configuration imagin- 
able, a one dimensional semi-infinite medium. The up- 
per picture in Fig. 2 shows the left end of an acoustic half 
space a; > with uniform dens ity pp and bulk modulus 
kq. The wave speed is cq = ^/kqJpo- Now replace the 
region < x < b with a shorter section 0<5 — a<a;<6 
filled with an acoustic metafiuid. The acoustic mirage ef- 
fect requires that an observer in a; > 5 hears a response as 
if the half space is as shown in the top of Fig. 2. This oc- 
curs if the metafluid region is such that: (i) no reflection 
occurs at the interface x = b, i.e. the acoustic impedance 
in the modified region is the same as before; and (ii) the 
round trip travel time of a wave incident from the right 
is unchanged. The impedance condition and the travel 
time requirement ensure the amplitude and phase of any 
signal is exactly the same as in the original half-space, 
and hence the mirage is accomplished. 

Let the acoustic metafluid have material properties p 
and K, with speed c = \fnjp. Conditions (i) and (ii) are 
satisfied if 



pc = poco, 



b 

Co 



(1) 



respectively, implying the density and bulk modulus in 
the shorter region are 



P = -Po, 
a 







(2) 



In this example the acoustic metafluid is another acoustic 
fluid, although with different density and bulk modulus. 
Note that the total mass of the metafluid region is un- 
changed from the original: po^ = pa- 



B. 2D mirage 

Consider the same problem in 2D under oblique wave 
incidence, Fig. 3. A wave incident at angle from 
the normal has travel time 6/(coCOS^o) in the original 
layer. If the shortened region has wave speed c, then 
the modified travel time is a/(ccos6') where is de- 
fined by the SncU-Descartcs law, ^sinfl ~ —smOn. At 
the same time, the reflectivity of the modified layer is 
R = {Z - Zn)/{Z + Zq) where Z = pc/ cos 0. The 
impedance and travel time conditions arc now 



pc 
cos 6 



PqCq 
cos 9q ' 



c cos 9 c cos ( 



(3) 



Solving for the modifled parameters implies n is again 
given by eq. (2)2 but the density is now 



Po 



a 2 cos^ 9q 



1 



l__sin2 2e„ 



(4) 



The mirage works only for a single direction of incidence, 
^0, and is therefore unsatisfactory. The underlying prob- 
lem here is that three conditions need to be met: Snell- 
Descartcs' law, matched impedances and equal travel 
times, with only two parameters, p and k. Some ad- 
ditional degree of freedom is required. 
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2. Solution of the 2D mirage problem 
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FIG. 3. The same as Fig. 2 but now for oblique incidence. 
The same strategy used for the ID case is no longer adequate; 
however a solution may be found if the metafluid in the layer 
of thickness a is allowed to have anisotropic inertia defined by 
the inertias p\ and pi in the x\ and directions. 



1. Anisotropic inertia 



The travel time is a/w • ei, and the impedance is now 
Z = piv/ cosO, where, referring to Fig. 3 ni = cos9 and 
122 = sin 9. Hence, the conditions for zero reflectivity and 
equal travel times are 



piv 



Poco 



cos 9 cos 9q 



Co cos 9o 



(9) 



Dividing the latter two relations implies k is given by 
(2)2. SncU-Descartes' law, v~^sm9 = Cq 
bined with eq. (9)i yields 



sm( 



com- 



pl 



(1-4) sin' 00, 



(10) 



while SncU-Dcscartes' law together with cq. (8)1 implies 

(11) 



cIpi 



fl-^Um'^o. 



P2' 



Comparison of (10) and (11) implies two identities for 
Pl and p2- In summary, the three parameters of the 
modified region are 



a a 

K = -rno, Pl = -Po, P2 = tPo- 
b a 



(12) 



These give the desired result: no reflection, and the same 
travel time for any angle of incidence. The metafluid 
layer faithfully mimics the wave properties of the original 
layer as observed from exterior vantage points. 



One method to resolve this problem is to introduce 
the notion of anisotropic mass density, see Fig. 3. The 
density of the metafluid region is no longer a scalar, but 
becomes a tensor: p ^ p. The equation of motion and 
the constitutive relation in the metafluid are 

pv = — Vp, p = — kV • V, (5) 

where v is particle velocity and p the acoustic pres- 
sure. Assuming two-dimensional dependence with con- 
stant anisotropic density of the form 



and eliminating v implies that the pressure satisfies a 
scalar wave equation 

P-CiP,ll-C2P,22 =0, (7) 

where Cj = ^/Ji/pj, j ~ 1,2. Equations (5) and (7) 
are discussed in greater detail and generality in §V, but 
for the moment we cite two results necessary for find- 
ing the metafluid in Fig. 3: the phase speed in direction 
n = niBi + 77-262 is V, and the associated wave or group 
velocity vector is w, where 

v'^ = c\n\ + 02^2, = v^^ [c^fiiei + £271262). (8) 



3. Anisotropic stiffness 

An alternative solution to the quandary raised by eq. 
(4) is to keep the density isotropic but to relax the 
standard isotropic constitutive relation between stress 
(T and strain e = [Vu -I- (Vu)*]/2 to allow for mate- 
rial anisotropy. Thus, the standard relation cr = —pl 
with p = K tr £ is replaced by the stress-strain relation 
for pentamode materials^'^ 

a = k{Q : e) Q, divQ = 0. (13) 

The physical meaning of the symmetric second order ten- 
sor Q is discussed later within the context of a more gen- 
eral constitutive theory. The requirement div Q = was 
flrst noted by Norris® and is discussed in §V. Standard 
acoustics corresponds to Q = I. 

Rewriting (13)i as cr = — pQ and using the divergence 
free property of Q, the equation of motion and the con- 
stitutive relation in the metafluid arc now 

pv = -QVp, p = kQ : Vv. (14) 

Eliminating v yields the scalar wave equation 

p- kQ : V(/9"^QVp) = 0. (15) 

General properties of this equation have been discussed 
by Norris^ and will be examined later in §V. For the 
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purpose of the problem in Fig. 3 Q is assumed constant 
of the form 



Q 



Qi 0' 
Q2 



(16) 



then it follows that the phase speed and group velocity 
in direction n are 



w ^v-\Clniei + C^n2e2), (17) 



where Cj = Qjy/n/p, j = 1,2. Proceeding as before, 
using Snell-Descartes' law and the conditions of equal 
travel time and matched impedance, yields 



a kq 
b K ' 



Q1Q2 



Since the important physical quantity is the product of k 
with Q (S) Q, any one of the three parameters k, Qi and 
Q2 may be independently selected. A natural choice is to 
impose div Q = at the interface, which means that the 
"flux" vector Qn is continuous, where n is the interface 
normal. In this case n = ei so that Qi = 1 and 



P 



-Po, 



-Kq, Q 



1 
A 



(18) 



Comparing the alternative metafluids defined by eqs. 
(12) and (18) note that in each case the density and the 
stiffness associated with the normal ei direction both 
equal their ID values given by eq. (2). The first metafluid 
of eq. (12) has a smaller inertia in the transverse direc- 
tion 62 (p2 < Pi)- The second metafluid deflned by eq. 
(18) has increased stiffness in the transverse direction 
{Q2 > Qi)- The net effect in each case is an increased 
phase speed in the transverse direction as compared with 
the normal direction (c2 > ci, C2 > Ci). 

The ID and 2D mirage examples illustrate the gen- 
eral idea of acoustic metafluids as fluids that replicate 
the wave properties of a transformed region. However, 
the methods used to flnd the metafluids are not easily 
generalized to arbitrary regions. How does one flnd the 
metafluid that can, for instance, mimic a full spherical 
region by a smaller shell? This is the cloaking problem. 
The key to the generalized procedure are the related no- 
tions of transformation and finite deformation, which are 
introduced next. 



III. THE TRANSFORMATION METHOD 

A. Preliminaries 

Let n and w denote the original and the deformed do- 
mains (the regions < X < b and 6 — a < a; < 6 in 
the examples of §11). The coordinates in each configu- 
ration arc X and x, respectively; the divergence opera- 
tors are Div and div, and the gradient (nabla) operators 
are V and V. Upper and lower case indices indicate 
components, Xj, Xi and the component form of div A is 
dAi/dxi = Ai^i or dAij/dxi = A^i when A is a vec- 
tor and a second order tensor-like quantity, respectively. 



and repeated indices imply summation (case sensitive!). 
Similarly DivB = Bijj. 

The finite deformation or transformation is defined by 
the mapping — > w according to x = x(X). In the 
terminology of finite elasticity X describes a particle po- 
sition in the Lagrangian or undeformed configuration, 
and X is particle location in the Eulerian or deformed 
physical state. The transformation x is assumed to be 
one-to-one and invertible^^. The deformation gradient is 
defined F = Vx with inverse F^^ = VX, or in com- 
ponent form Fii ~ dxi/dXj, Ff^^ ~ dXj/dxi. The 
Jacobian of the finite deformation is A = det F ~ |F|, 
or in terms of volume elements in the two configura- 
tions, A = dv/dV. The polar decomposition implies 
F = VR, where the rotation R is proper orthogonal 
(RR* = R'R = I, detR = 1) and the left stretch tensor 
V gSym+ is the positive definite solution of = FF*. 
Note for later use the kinematic identities'^ 

div (A^'F) = 0, Div (AF^^) = 0, (19) 

and the expression for the Laplacian in X in terms of 
derivatives in x, i.e., the chain rule® 



V / = A div (A-^VV/) 



(20) 



B. The transformation method 



The undeformed domain Q is of arbitrary shape and 
comprises a homogeneous acoustic fluid with density pQ 
and bulk modulus kq. The goal is to mimic the scalar 
wave equation in 



(ko/po)V p = 0, xen, 



(21) 



by the wave equation of a metafluid occupying the de- 
formed region uj. The basic result^'*'-'^^ is that eq. (21) is 
exactly replicated in lu by the equation 

p- Kdiv (p^'Vp) = 0, xGuj, (22) 

where the bulk modulus and inertia tensor are 



K = KqA, 



PoAV- 



(23) 



The equivalence of eq. (21) with eqs. (22) and (23) 
is evident from the differential equality (20). The idea 
is to use the change of variables from X to x to iden- 
tify the metafluid properties. Equation (23) describes a 
metafluid with anisotropic inertia and isotropic elastic- 
ity. It can be used for modeling acoustic cloaks but has 
the unavoidable feature that the total effective mass of 
the cloak is inflnite. This problem, discussed by Norris®, 
arises from the singular nature of the flnite deformation 
in a cloak which makes AV~^ non-integrable. This type 
of fluid, which could be called an inertial fluid, appears 
to be the main candidate considered for acoustic cloak- 
ing to date. The major exception is Milton et al.^' who 
considered fluids with properties of pentamode materials, 
although as we will discuss in Section IV, their findings 
are of limited use for acoustic cloaking. 
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1. Pentamode materials 

Norris® showed that eq. (20) is a special case of a more 
general identity: 



V^p = AQ : V{A-^q-^V^S/p), 



(24) 



where Q is any symmetric, invertible and divergence free 
(div Q = 0) second order tensor. The increased degrees of 
freedom afforded by the arbitrary nature of Q means that 
(21) is equivalent to the generalized scalar wave equation 



m w, 



P 



kQ : V(p"^QVp) 0, xGw, 



(25) 



where the modulus k and the inertia follow from a com- 
parison of eqs. (21), (24) and (25) as 



kqA, 



p = poAQ VQ. 



(26) 



As will become apparent later, these metafluid parame- 
ters describe a pentamode material with anisotropic in- 
ertia. For the moment we return to the acoustic mirages 
in light of the general transformation method. 



C. Mirages revisited 

The 2D mirage problem corresponds to the following 
finite deformation xi = b — a + ab^^Xi, X2 = X2 for 
< ATi < 6, — 00 < X2 < 00. The deformation gradient 
is 



f 
1 



(27) 



implying R = I,V = 
therefore implies 



K = -Ko, 



F and A = a/b. Equation (23) 



>o 
tPo 



(28) 



These agree with the parameters found in §11, eq. (12). 

Using the more general formulation of (25) and (26) 
with the arbitrary tensor chosen as Q = A^^V yields the 
metafluid described by eq. (18). It is interesting to note 
that although p of eq. (26) is, in general, anisotropic, it 
can be made isotropic in this instance by any Q propor- 
tional to V. Keeping in mind the requirement seen above 
that Qi = 0, we consider as a second example of (25) the 
case Q = A^^V. The mirage can then be achieved with 
material properties 



P 



-Po, 



-Kq, Q 



1 
A 



(29) 



This again corresponds to a pentamode material with 
isotropic inertia, equal to that of (18). These two ex- 
amples illustrate the power associated with the arbitrary 
nature of the divergence free tensor Q. There appears to 
be a multiplicative degree of freedom associated with Q 
that is absent using anisotropic inertia. As will be ev- 
ident later, this degree of freedom is related to a gauge 
transformation. 



IV. THE iVIOST GENERAL TYPE OF ACOUSTIC 
IVIETAFLUID 

A. Summary of main result 

In order to make it easier for the reader to assimilate, 
the paper's main result is first presented in the form of a 
theorem. In the following Qq and Q are arbitrary sym- 
metric, invertible and divergence free second order ten- 
sors. 

Theorem 1 The kinetic and strain energy densities in 
O, of the form 

To = U*PoU, Wo = Ko{Qo : VU)', (30) 
respectively, are equivalent to the current energy densities 



where 



r = u*pu, ly = «;(Q:Vu)^ (31) 

K = Ako, (32a) 
p = AQF-*Qo-VoQo 'F-iQ. (32b) 



Discussion of the implications are given following the 
proof. 



B. Gauge transformation 

The energy functions per unit volume in the unde- 
formed configuration, Tq and Wq, depend upon the in- 
finitesimal displacement U in that configuration. The 
kinetic energy is defined by the density Pq, while the 
strain energy is Wo = CoijklUjjUl,k, where Cqijkl 
are elements of the stiffness. The density and stiffness 
possess the symmetries poij = poji, Cqijkl = Coklij, 
Cqijkl = Cojikl- The total energy is Eo = To+Wq and 
the total energy E — T + W per unit deformed volume 
is, using EodV ~ E dv, simply E = A^^Eo- 

Our objective is to find a general class of material pa- 
rameters {pq, Co} that maintain the structure of the en- 
ergy under a general transformation x- Structure here 
means that the energy remains quadratic in velocity and 
strain. In order to achieve the most general form for the 
transformed energy introduce a gauge transformation for 
the displacement. Let 



U 



G*u, 



(33) 



or Ui = UiGii in components, where G is independent of 
time but can be spatially varying. Thus, using the chain 
rule Ujj ~ UjjFii, yields E ~ T + W, where the kinetic 
and strain energy densities are 



T = u*pu. 



(34a) 



W - A-^CijKLFuFkK{Gjju,)4GiLUi)^k, (34b) 
and the transformed inertia is 



p = A-'GpoGK 



(35) 
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The kinetic energy has the required structure, quadratic 
in the velocity; the strain energy however is not in the 
desired form. The objective is to obtain a strain energy 
of standard form 

W = Ci-jkiUj.m.k, (36) 

where C has the usual symmetries: Cijki = CkHj, Cijki = 
Cjiki ■ 

The question of how to convert W of (34b) into the 
form of (36) for generally anisotropic elasticity stiffness 
Co will be discussed in a separate paper. The objective 
here is to find the largest class of metafluids that includes 
all those previously found. 

C. Pentamode to pentamode 

In order to proceed assume that the initial stiffness 
tensor is of pentamode form^'^ 

Co = AtoQo ® Qo: Div Qo = 0, (37) 

that is, CoijKL = KoQoijQoKL where Qq = Qo is a pos- 
itive definite symmetric second order tensor. The tensor 
Qo is necessarily divergence frce^. 

The strain energy density in the physical space after 
the general deformation and gauge transformation is now 

W^KoA-'[Qoij{u,G,j)j]\ (38) 

Consider 

Qoij{u]Gjj),i = QoijGjjUjj + QoijGjjjUj. (39) 

In order to achieve the quadratic structure of (36) the 
final term in (39) must vanish. Since u is considered 
arbitrary this in turn implies that QoijGjjj must vanish 
for all j. With no loss in generality let 

G* = Qo^P, (40) 

or Gjj = QojM-P^^j components. Then using the iden- 
tity for the derivative of a second order tensor. (A^^).q, = 
— A^^A qA^^ where a can be any parameter, gives 

= QoijGjjj 
= -QoijQqjkQokljQoImPmj + PljJ 
= -QoiljQqImPmj + Plj,I 
- PijJ, (41) 

where the important property (37)2 has been used. 
Hence, Div P = 0. Then using (39) yields 

QoijiujGjj)j = PijUjj 

— Plj Pi I^j^i 

= AQyUj-i, (42) 
where the tensor Q is defined by 

P = AF"iQ ^ Q = A-iFP. (43) 



The condition DivP = becomes, using (19)2, 

Pi J J = {AFj-^'Q,,)j = AF^^Q^jj = AQ„-,„ (44) 
implying 

divQ = 0. (45) 

It has therefore been demonstrated that if the gauge 
transformation is defined by 

G*=AQo1f-1Q, (46) 

where Q is divergence free in physical coordinates, then 
the strain energy (38) is = KoA(^QijUj_i) . This is of 
the desired form. eq. (36), with Cijki = i^-QijQki, hence 
completing the proof of Theorem 1 . 

D. Discussion 

1. Equivalence of physical quantities 

Theorem 1 states that the pentamode material defined 
by stiffness kq with stress-like tensor Qg and anisotropic 
inertia is converted into another pentamode mate- 
rial with anisotropic inertia. The properties of the new 
metafluid are defined by the original metafluid and the 
deformation-gauge pair {F, G} where F is arbitrary and 
possibly inhomogeneous, and G is given by eq. (46) with 
Q symmetric, positive definite and divergence free but 
otherwise completely arbitrary. The special case of a fluid 
with isotropic stiffness but anisotropic inertia, eq. (23), 
is recovered from Theorem 1 when the starting medium 
is a standard acoustic fluid and Q is taken to be I. 

It is instructive to examine how physical quantities 
transform: we consider displacement, momentum and 
pscudoprcssurc. Eliminating G, it is possible to express 
the new displacement vector in terms of the original, 

u = A-iQ-^FQgU. (47) 

Physically, this means that as the metafluid in uj acous- 
tically replicates that in CI, particle motion in the latter 
is converted into the mimicked motion defined by (47). 

Define the momentum vectors in the two configura- 
tions, 

mo = PqU, m = pu. (48) 

Equations (35) and (46) imply that they are related by 

m = QF-*Qoimo. (49) 

The transformation of momentum is similar to that for 
displacement, eq. (47), but with the inverse tensor, i.e. 
u = G"*U while m = A^^Gmo. 

Stress in the two conflgurations is defined by Hooke's 
law in each: 

o-o = Co : VU, cr = C : Vu, (50) 
where Co and C arc fourth order elasticity tensors, 

Co = KoQo «) Qo, C = KQ(g)Q, (51) 
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that is, Cijki = nQijQki, etc. Using (37), (42) and (51), 
yields 

o-o = -pQo, cr = -pQ, (52) 

where p is the same in each configuration, 

p = -kQ : Vu = -KoQo : VU. (53) 

The quantity p is similar to pressure, and can be exactly 
identified as such when Q is diagonal, but it is not pres- 
sure in the usual meaning. For this reason it is called the 
pseudopressure. It is interesting to compare the equal 
values of p in r2 and uj with the more complicated rela- 
tions (47) and (49) for the displacement and momentum. 

2. Equations of motion 

The equations of motion can be derived as the Euler- 
Lagrange equation for the Lagrangian T — W^. A succinct 
form is as follows, in terms of the the momentum density 
m and the stress tensor cr: 

m = Vcr. (54) 

The constitutive relation may be expressed as an equa- 
tion for the pseudopressure , 

P=-kQ:Vu, (55) 

while eqs. (53), (48) and (54) imply that the acceleration 
is 

ti = -p^Q\/p. (56) 

Eliminating p between the last two equations implies that 
the displacement satisfies 

kQQ : Vii - pii = 0. (57) 

This is, as expected, the elastodynamic equation for a 
pentamode material with anisotropic inertia. Alterna- 
tively, eliminating ii between eqs. (55) and (56) yields a 
scalar wave equation for the pseudopressure, 

p- kQ : V(p"iQVp) = 0. (58) 

This clearly reduces to the standard acoustic wave equa- 
tion when Q = I and p is isotropic. 



Using (19)i and (37)2 yields Qy,* = ^'^QoijFjjj. 
Hence, this particular Q can only satisfy the requirement 
(45) that div Q = if 

Milton et al.^^ considered Qq isotropic (diagonal), in 
which case (60) means that the only permissible finite de- 

2 

formations are harmonic, i.e. those for which V x = 0. 
In short, the negative findings of Milton et al.^"'^ are a 
consequence of constraining the gauge to G = F, which 
in turn severely restricts the realizability of metafluids 
except under very limited types of transformation defor- 
mation. The main difference in the present analysis is 
the inclusion of the general gauge transformation which 
enables us to find metafluids under arbitrary deforma- 
tion. 

V. PROPERTIES OF ACOUSTIC METAFLUIDS 

The primary result of the paper, summarized in The- 
orem 1, states that the class of acoustic metafluids is 
defined by the most general type of pentamode mate- 
rial with elastic stiffness kQ (E) Q where divQ = 0, and 
anisotropic inertia p. We now examine some of the un- 
usual physical properties, dynamic and static, to be ex- 
pected in acoustic metafluids. Some of the dynamic prop- 
erties have been discussed by Norris®, but apart from 
Milton and Cherkaev^" no discussion of static effects has 
been given. 



A. Dynamic properties: plane waves 

Consider plane wave solutions for displacement of the 
form u(x, i) = qe'^-'C^'^-"*), for |n| = 1 and constant q, 
k and v, and uniform metafluid properties. Non triv- 
ial solutions satisfying the equations of motion eq. (57) 
require 

(K(Qn)®(Qn)-pt;2)q = 0. (61) 

The acoustical or ChristoffcP° tensor K(Qn) (g) (Qn) is 
rank one and it follows that of the three possible solu- 
tions for w^, only one is not zero, the quasi-longitudinal 
solution 



3. Relation to ttie findings of Milton et al 

The present findings appear to contradict those of Mil- 
ton et al.^^ who found the negative result that it is not 
in general possible to find a metafluid that replicates a 
standard acoustic medium after arbitrary finite deforma- 
tion. However, their result is based on the assumption 
that G = F (their eq. (2.2)). Equation (46) implies that 
Q must then be 

Q = A-iFQoF*. (59) 



w^ = Kn Qp ^Qn, q = p ^Qn. (62) 

The slowness surface is an ellipsoid. The energy fiux 
velocity^" (or wave velocity or group velocity ) is 

w = u^^KQp^iQn. (63) 

v^r is in the direction Qq, and satisfies w • n = u, a 
well known relation for generally anisotropic solids with 
isotropic density. 
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B. Static properties 

1. Five easy modes 

The static properties of acoustic metafluids are just as 
interesting, if not more so. Hooke's law (52)2 is 

cr = Cs, (64) 

where e = e* is strain and the stiffness C is defined by 
eq. (50)2- The strain energy is = k(Q : e)^. Note 
that C is not invertible in the usual sense of fourth or- 
der elasticity tensors. If considered as a 6 x 6 matrix 
mapping strain to stress then the stiffness is rank one: it 
has only one non-zero eigenvalue. This means that there 
are five independent strains each of which will produce 
zero stress and zero strain energy, hence the name pen- 
tamode^^ . The five "easy" strains are easily identified in 
terms of the principal directions and eigenvalues of Q. 
Let 

Q = AiQiqi + A2q2q2 + Asqgqg, (65) 

where {qj^, qj, q3} is an orthonormal triad. Three of the 
easy strains are pure shear: q^q^ + q^q,, i 7^ j and the 
other two are A3q2q2—A2q3q3 and Aiq2q2 — A2qiqi. Any 
other zero-energy strain is a linear combination of these. 
Note that there is no relation analogous to (64) for strain 
in terms of stress because only the single "component" 
Q : £ is relevant, i.e., energetic. 

It is possible to write cr in the form (52) where p = 
— kQ : £. Under static load in the absence of body force 
choose K such that p =constant, or equivalently, eq.(37)2. 
The relevant strain component is then Q : e = —K~^p 
and the surface tractions supporting the body in equi- 
librium are t = a ■ n = — pQ • n. Figure 4 illustrates 
the tractions required to maintain a block of metafluid 
in static equilibrium. Note that the traction vectors act 
obliquely to the surface, implying that shear forces are 
necessary. Furthermore, the tractions are not of uniform 
magnitude. These properties are to be compared with a 
normal acoustic fluid which can be maintained in static 
equilibrium by constant hydrostatic pressure. 



2. div Q = 

The gedanken experiment of Fig. 4 also implies that Q 
has to be divergence free. Thus, imagine a smoothly vary- 
ing but inhomogeneous metafluid C = kQ g) Q in static 
equilibrium under the traction t = —pQ ■ n for constant 
p. The divergence theorem then implies div Q = every- 
where in the interior. This argument is a bit simplistic 
- but it provides the basis for a more rigorous proof^. 
Thus, stress in the metafluid must be of the form — /Q 
where Q is a scalar multiple of Q. Local equilibrium re- 
quires div/Q = 0, or Vln/ = — Q divQ. This can 
be integrated to find / to within a constant. Now de- 
fine Q = /Q, and note that the tractions must be of the 
form t = pQ ■ n for constant p. The normalized Q is 
divergence free. 




FIG. 4. A rectangular block of metafluid is in static equilib- 
rium under the action of surface tractions as shown. The two 
orthogonal arrows inside the rectangle indicate the principal 
directions of Q (30° from horizontal and vertical) and the 
relative magnitude of its eigenvalues (2:1). The equispaced 
arrows faithfully represent the surface loads. 



3. Non-horizontal free surface 

Consider the same metafluid in equilibrium under a 
body force, e.g. gravity. Assuming the inertia is isotropic 
(cf. the comments about inertia at zero frequency in 
§VI), 

div cr + pg = 0. (66) 

Use eq. (52)2 with divQ = and the invertibility of Q, 
implies 

Vp = pQ-ig. (67) 

For constant Q~^/Og this can be integrated to give an 
explicit form for the pseudopressure, 

p=(x-xo)-Q-Vg, (68) 

where xq is any point lying on the surface of zero pres- 
sure. Unlike normal fluids, the surface where p = does 
not have to be horizontal, see Fig. (5). The pseudopres- 
sure increases in the direction of g, as in normal flu- 
ids. However, it is possible that p varies in the plane 
X • g = 0. For instance, the traction along the lower 
surface in Fig. (5) decreases in magnitude from left to 
right. 

VI. IVIETAFLUIDS WITH ISOTROPIC DENSITY 
A. Necessary constraints on the finite deformation 

The most practical case of interest is of course where 
the initial properties are those of a standard acoustic fluid 
with isotropic density and isotropic stress. The circum- 
stances under which the mapped inertia p is also isotropic 
are now investigated. Acoustic metafluids with isotropic 
inertia are an important subset since it can be argued 
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FIG. 5. The same metafluid of Fig. 4 with isotropic density 
is in equilibrium under gravity. The upper surface is traction 
free but non-horizontal, an essential feature of metafluids. For 
this particular metafluid the top surface makes an angle of 
19.11° with the horizontal. Also, the tractions on the hori- 
zontal lower boundary are inhomogeneous although parallel. 
The small rectangle is discussed in Fig. 6. 



that achieving anisotropic inertia could be more difficult 
than the anisotropic elasticity. Indeed, the very concept 
of anisotropic inertia is meaningless at zero frequency, 
unlike anisotropic stiffness. 

Assuming Pq = p^l and Qq = I then the current den- 
sity becomes, using eq. (32b) and the fact that Q is 
symmetric, 



p = poAQV 2q. 



(69) 



If p = pi then eq. (69) implies Q must be of the form 
Q = pi/2(poA)-i/2v. (70) 

Thus, Q is proportional to the stretch tensor V and the 
coefficient of proportionality defines the current density. 
It is not in general possible to choose Q in the form 



Q = aV, 



(71) 



where a ^ 0. Certainly, Q of (71) is symmetric and 
invertible but not necessarily divergence free. The lat- 
ter condition requires div(aV) = which in turn may 
be expressed Vina = — V~^divV. The necessary and 
sufficient condition that V~^divV is the gradient of a 
scalar function, and hence a can be found which makes 
Q of (71) possible, is that V satisfy 




FIG. 6. This shows a possible microstructure for a micro- 
scopic rectangular region of the metafluid in Figs. 4 and 5, 
see the latter. The sides of the rectangle are aligned with the 
principal axes of Q. The metafluid is a pentamode material 
comprising a regular array of small beads such that each is in 
lubricated point contact with its three neighbors (2D). The 
dashed lines indicate the directions of the forces acting be- 
tween the small oval-shaped beads. Although the structure 
as shown is unstable under shear, a realistic metafluid might 
contain some stabilizing mechanisms to enhance its rigidity. 
Details will be provided in a forthcoming paper. 



This condition is not very useful. It does, however, indi- 
cate that the possibility of achieving isotropic p depends 
on the underlying finite deformation; there is a subset 
of general deformations that can yield isotropic inertia. 
The deformation gradient F has nine independent ele- 
ments, V has six, and R has three. The condition (72) is 
therefore a differential constraint on six parameters. We 
now demonstrate an alternative statement of the condi- 
tion div(aV) = in terms of the rotation R. This will 
turn out to be more useful, leading to general forms of 
potential deformation gradients. 

Substituting F ~ a^^QR into the identity (19)i and 
using (45) implies that 



(73) 



Using Q = aFR* and the identity F*V = V along with 
a ^ 0, yields 



RjM{A-^a-^RjK) 



0. 



Then using the identity (RjKRjAi) j^j 
yields 



0, eq. 



/9, 



K 



RjkRj 



M.M 



R*DivR*, 



(74) 
(74) 

(75) 



where /3 = — In(Aa). 

The necessary and sufficient condition that (75) can be 
integrated to find /3 is V A V/3 = 0, or using (75) 



CurlR*DivR* = 0. 



(76) 



curlV^MivV = 0. 



(72) 



The integrability condition (76) is in general not satisfied 
by R, except in trivial cases. Norris® noted that isotropic 
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density can be obtained if R = I. This corresponds to 
(3 =constant, and it can be realized more generally if R 
is constant. Hence, 

Lemma 1 // the rotation R is constant then a normal 
acoustic fluid can be mapped to a unique metafluid with 
isotropic inertia: 

= Ko(DivU)VpoU-U inn, (77) 

is equivalent to the current energy density 

E ^ \{Y -.^xif + ph-ii into, (78) 

where 

X = A-^Ko, p = A-Vo. (79) 

The total mass of the deformed region to is the same as 
the total mass contained in n. 



Lemma 2 // the rotation R satisfies (76), or equiva- 
lently, if 9 and a. satisfy the condition (82), then a normal 
acoustic fluid can be mapped to a unique metafluid with 
isotropic inertia according to eqs. (77) and (78) with 

A = A-^e-^^Ko, p = K^^e-'^^p(,, (83) 

where the function (3 is defined by eq. (75). 



C. Simplification in 2D 

The integrability condition (76) simplifies for the im- 
portant general configuration of two-dimensional spatial 
dependence. In this case a is constant, 6 = ^(Xj^) where 
Xj^ • a = 0. Equation (82) then reduces to 

V^e ^ 0, (84) 



The parameter A is used to distinguish it from hl = Akq, 
because in this case Q = A^^V. Also, the displacement 
fields are related simply by u = RU. 

As an example of a deformation satisfying Lemma 1: 
X = /(X • AX) AX for any constant positive definite 
symmetric A. This type of finite deformation includes 
the important cases of radially symmetric cloaks. Thus, 
Norris^ showed that radially symmetric cloaks can be 
achieved using pcntamodc materials with isotropic iner- 
tia. 



B. General condition on the rotation 

The results so far indicate that isotropic inertia is 
achievable for transformation deformations with constant 
rotation. We would however like to understand the 
broader implications of eq. (76). The rotation can be 
expressed in Euler form 



R = exp {6 axt a) 



(80) 



where 9 is the angle of rotation, the unit vector a is the 
rotation axis, and the axial tensor axt(a) is a skew sym- 
metric tensor defined by axt(a)b = a A b. The vector 
Oa. encapsulates the three independent parameters in R. 
The integrability condition (76) is now replaced with a 
more explicit one in terms of ^?(X) and a(X). It is shown 
in the Appendix that 



R* DivR* 



aA Ve- 



(81) 



1. Example 

Consider finite deformations with inhomogeneous ro- 
tation 



= 6*0+7X1, a = e3, 



(85) 



for constants 6*0 and 7. This satisfies (84) and there- 
fore eq. (75) can be integrated. The metafiuid in uj 
has isotropic density and pentamode stiffness given by 
Lemma 2, where (3 — 7(^2 — X02)). The constants 6*0 
and Xq2 may be set to zero, with no loss in generality. 

As an example of a deformation that has rotation of 
the form (85), consider deformation of the region Vl = 
[— ^] X [0, L] X K according to 







All 


Ar2 


0' 






^12 


A22 















a 



7^1 sin 7X1 e^^^'^ 



cos 7X1 e"'*'-^") 

^3 



(86) 

where a > and the 2x2 symmetric matrix A with ele- 
ments Aij is positive definite. The deformation gradient 
is F = V(X2)R(Xi) with V = Ae''^^^ + aesea, and 



R 



cos 7X1 — sin 7X1 
sin 7X1 cos 7X1 
1 



(87) 



The mapped metafluid is defined by the energy density 
in u} 



where the vector Z follows from eq. (A6). In particular, 
Z vanishes if the axis of rotation a is constant. In general, 
for arbitrary spatial dependence, (81) implies that the 
integrability condition (76) is equivalent to the following 
constraint on the rotation parameters: 



aV'e- 



-(a-V)V6l+(Ve'-V)a-(V-a)V6'+CurlZ 
In summary. 



0. (82) 



E 



1 



a det A 



Ko[(A : Vu) 



2 -27A-2 



+ (au3,3)^] +POU-U 



In particular, it has isotropic inertia. 

This example is not directly applicable to modeling a 
complete acoustic cloak. However, it opens up the possi- 
bility of patching together metafluids with different local 
properties, each with isotropic inertia so that the entire 
cloak has isotropic mass density. 
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VII. SUMMARY AND CONCLUSION 

Whether it is the simple ID acoustic mirage of Fig. 2 
or a three-dimensional acoustic cloak, we have seen how 
acoustic stealth can be achieved using the concept of do- 
main transformation. The fluid in the transformed re- 
gion exactly replicates the acoustical properties of the 
original domain. The most general class of material that 
describes both the mimic and the mimicked fluids is de- 
fined as an acoustic metafluid. A general procedure for 
mapping/ transforming one acoustic metafluid to another 
has been described in this paper. 

The results, particularly Theorem 1 in §1V, show that 
acoustic metafluids are characterized by as few as two 
parameters (p, ^/k) and as many as twelve (p, y^Q). 
This broad class of materials can be described as pcn- 
tamode materials with anisotropic inertia. It includes 
the restricted set of fluids with anisotropic inertia and 
isotropic stress (Q = I). 

The arbitrary nature of the divergence free tensor Q 
adds an enormous amount of latitude to the stealth prob- 
lem. It may be selected in some circumstances to guar- 
antee isotropic inertia in cloaking materials, examples of 
which are given elsewhere*. In this paper we have de- 
rived and described the most general conditions required 
for p to be isotropic. The conditions have been phrased 
in terms of the rotation part of the deformation, R. If 
this is a constant then the cloaking metafluid is defined 
by Lemma 1. Otherwise the condition is eq. (76) with the 
metafluid given by Lemma 2. The importance of being 
able to use metafluids with isotropic inertia should not 
be underestimated. Apart from the fact that it resolves 
questions of infinite total effective mass* isotropic iner- 
tia removes frequency bandwidth issues that would be 
an intrinsic drawback in materials based on anisotropic 
inertia. 

This paper also describes, for the first time, some of the 
unusual physical features of acoustic metafluids. Strange 
effects are to be expected in static equilibrium, as illus- 
trated in Figs. 4 and 5. These properties can be best 
understood through realization of acoustic metafluids, 
and a first step in that direction is provided by the type 
of microstructure depicted in Fig. 6. The macroscopic 
homogenized equations governing the microstructure are 
assumed in this paper to be those of normal elasticity. 
It is also possible that the large contrasts required in 
acoustic metafluids could be modeled with more sophis- 
ticated constitutive theories, such as non-local models or 
theories involving higher order gradients. There is con- 
siderable progress to be made in the modleing, design 
and ultimate fabrication of acoustic metafluids. 

In addition to the degrees of freedom associated with 
the tensor Q, the properties of metafluids depend upon 
the finite deformation through V. Even in the simple ex- 
ample of the ID mirage, one could arrive at the lower pic- 
ture in Fig. 2 through different finite deformations. This 
raises the question of how to best choose the nonunique 
deformation gradient F. The present results indicate 
some strategies for choosing F to ensure the cloak in- 
ertia has isotropic mass, and the cloaking properties are 
in effect determined by the elastic pentamodal material. 



Li and Pendry^ consider other optimal choices for the fi- 
nite deformation. Combined with the enormous freedom 
afforded by the arbitrary nature of Q, there are clearly 
many optimization strategies to be considered. 
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APPENDIX A: DERIVATION OF EQ. (81) 

Equation (80) can be written in Euler-Rodrigues 
form^^ 



R = I + sin6'S + (l-cos6')S^ 



where 



axt(a) = 



-as 02 
03 — fli 
—02 fli 



(Al) 



(A2) 



Equation (Al) can be used, for instance, to find 9 = 
(1 — trR)] and hence a from S = (R — 



cos 



-Irl 



R )/(2sin0). For the sake of notational brevity, in the 
remainder of the Appendix c = cos^ and s = sin 9. 
Explicit differentiation of (Al) yields 

R^JJ = {cS^J - sP,j)9j + sSuj - (1 - c)Pu.j, (A3) 

where P = -S^ = I - aa. Noting that Sijj = -(V A a)i 
and Pij.j = — ajV • a — a • Va^, implies 



DivR' = ca A V6' - sV A a - s(I - aa) • \79 
+ (l-c) [a(V • a) + (a • V)a] . 



(A4) 



Multiplying by R* using (Al), gives after some elimina- 
tion and simplification, 



R*DivR* 



aLAV9 + Z, 



(A5) 



where 



Z = -s[cl + (1 - c)aa] • V A a (1 - c)a(V • a) 

- (1 - c)[I- (1 - c)aa] • (a- V)a 

- s(l - c)aA (a- V)a. (A6) 
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